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Abstract The study area was located in I inneshin Natural Rcscrxe. Xaoxing’an Mountains. Northeastern China. Korean 
pine forests are the Apical forest cross stems and landscapes in this region. It is a high degrees of spatial and temporal het¬ 
erogeneity at different settles, which effected on landscape pattern and processes. In this paper we used the data of 144 plots 
and semis ariogrnm It' analyze spatial heterogeneity of old growth forests of Korean pine in landscape lesel. Model for forest 
sarialions by isotropic semisariogram is linear with sill flic spatial heterogencits is dependent on settles and directions in 
Korean pine forests. Patterns of forest types in space were resulted from complex interactions between physical and biologi¬ 
cal forces We used 20 metres for interpolation inters al to estimate (he sallies of unsampled area. Comparing the results with 
held data, block kriging and mapping are an effectts e techniques to simulate landscape pattern. 
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Introduction 

Landscape can be considered a spatiall\ heterogeneous 
area, and spatial pattern refers to the arrangement of 
structural heterogeneity of environments, resource and 
biological systems in space (Forman el til 1086. I timer 
1087. 1080. Kisser 1087. I timer el al 1001). Turner 
(1087. 1080) considered that spatial pattern has impor¬ 
tant consequences regarding the effective dispersal of 
organisms among landscape elements and the spread of 
disturbance across (he landscape. In landscape ecology, 
spatial patterns xxere caused In different natural proc¬ 
esses and interactions under differenl scales, and xxere 
distributed regularity in landscape mosaic. I 'here for. ii 
is one of the main questions in landscape ecological 
studies(I.i el al 1088. Wiens 1088. Turner 1080). There 
exists a serlieal and hori/onlal satiations of segeta- 
tionle.g.. lores! stands) and other resource (c.g.. soil) in 
landscape les el. I hese sarialions are communis in 
landscape and have txvo components: one is the struc¬ 
tural components.g.. patch shapes, si/e) xxhich can he 
identified, another is the random xxhich is unknown and 
uninterprelediBurrotigh 1085). Both ol’xxhich are the 
center in landscape studies) Forman el al 1086). 

l andscape slruclure must lie identified and quantilied 
in meaningful ways before the interactions between 


landscape pattern and ecological processes can be un¬ 
derstood!Turner 1080. |00|). There are mans quantita¬ 
tive wavs for spatial pattern analyzed (Routine 1082. 
O’Neill el al 1088. Turner 1088. 1001. Milne 1088. 
Burrough 1086. (iardner et al 1087) and spatial simula¬ 
tion of landscape pattern and changes (Turner 1087. 
1088. Bartrcll et al 1001. Merriam el al I00|. Gardner 
el al 1001). which were developed during p;ist over 
deeatle. Spaliai slalisties. e.g. spatial autocorrelation, is 
one of the spatial models. Semis ariogram. as spatial 
autocorrelation, examines spatial variation and correla¬ 
tion of natural phenomenal Mathcron 1065. .lourncl el al 
1078.Webster 1085. CTessie 1001). Semivariogram is 
based on the theory of regionalized variables and sshile 
recent applications of xxhich in ecological research base 
demonstrated its effectiveness to summarize spaliai 
data (Robertson 1087. Robertson el al 1088. Palmer 
1088. Fortin et al 1080. Legendre et al 1080. Levin et al 
1080. Rossi et al 1002). One of the prime reasons lor 
obtaining a semivariogram is to use For spatial estima¬ 
tion. e.g. kriging (Journel et al 1078. Webster 1085). 
Kriging is a technique of making optimal, unbiased 
estimates ol regionalized variables at unsampled loca¬ 
tions using the structural properties of the semivari¬ 
ogram and the initial set of data values)Journel et al 
1078. I rangmar et al 1085). Ihcrefor. kriging can be 
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used for landscape simulation and mapping spatial pat- 
tem(e.g., patches) by landscapefor spatial) data. 

Korean pine(Pinus koraiensis) forests, which often 
mixed with hardwoods, e.g. birch (Belula costata ), 
ash(Fraxinus mandshurica ), walnut (.Juglans man- 
shurica ), corktr ee(Phellodendron amuren.se), linden 
(Tilia amurensis ), typical old growth forest ecosys- 
tem(Chen 1982, Li et al. 1993, Wang et al. 1996), is a 
main important forest landscape in Northeast Parts of 
China. It is high degrees of spatial and temporal hetero¬ 
geneity at different scales(Wang et al 1997), which ef¬ 
fected on landscape pattern and processes. There are 
some different landscape patterns in Korean pine for¬ 
ests. The distribution of species, structures of age, en¬ 
ergy, nutrient cycling, and regeneration would be re¬ 
lated to the sizes, shapes of patterns. But, studies of 
landscape ecology, specifically, spatial pattern for Ko¬ 
rean pine forest have been lacking during the past. We 
should know the development and dynamics of spatial 
heterogeneity, interactions and exchanges across het¬ 
erogeneous landscape, the influences of spatial hetero¬ 
geneity on .biotic and a biotic processes, and the man¬ 
agement of spatial heterogeneity(Turner 1989) at dif¬ 
ferent scales in landscape of Korean pine forests, we 
may understand well the trues of structures, functions 
and dynamics of Korean pine forest ecosystems in 
landscape level. 

In this paper we use semivariogram model to analyze 
spatial heterogeneity and block kriging(Journel et al 
1978, Webster 1985, Trangmar et al 1985) to simulate 
landscape pattern of Korean pine forests, and discuss 
the relations between spatial heterogeneity and patterns. 
We specifically discuss (1) the characteristics, degrees 
and scales of spatial heterogeneity; (2)kriging and 
mapping for landscape patterns, and (3)the sizes and 
shapes of spatial patterns of Korean pine forests. 

Theory of Block Kriging 

I he method of estimation embodied in regionalized 
variable theory is known in earth sciences as kriging, 
after Krigef 1966), who first devised it empirically for 
use in the South African goldfields. It is essentially a 
means of weighted local averaging in which the 
weights are chosen so as to give unbiased estimates 
while at the same time minimizing the estimation vari¬ 
ance. Kriging is this sense optimalf Webster 1985). 

In block kriging, a value for an area or block with its 
center at A„ is estimated. The kriged value of property 
Z for any block V is a weighted average of the ob¬ 
served values X-, in the neighborhood of the block, i.e.. 


Z(/) = ^/ i Z(.V,) (1) 

/=! 

where n is the number of neighboring samples Z(XJ 
and K are weights applied to each Z(Xj. The weights 
are chosen so that the estimate z(V) of the true value 
Z(V) is unbiased, i.e., 

/;[Z(n-Z(L)] = 0 (2) 

and the estimation variance a J for block V is mini¬ 
mized, i.e., 

a" = VAR[Z(L)-Z(( .)] = minimum (3) 

The weights placed on each neighboring sample sum to 
1, and their unique combination for which c - is mini¬ 
mized can be obtained when 

li 

^A/yt.V,, n + p,. =y(L, n (4) 

/ I 

n 

and unbias conditions is YW = j. The estimated vari- 

i i 

ance is shown in below 

ol = n + fL.-YflA l ) (5) 

/- I 

where y(\ i d') is the average semivariance between 
the sample points X, in the neighborhood and those in 
the block V. y(F.F) is the average semivariance be¬ 
tween all points within V (i.e., the within block vari¬ 
ance of classical statistics)(Journel et al 1978), and // v 
is the Lagrangian parameter associated with the mini¬ 
mization. Eq.(4) is the n+1 equations of the kriging 
system. Kriging thus provides not only unbiased esti¬ 
mates of minimum variance, but also a measure Of the 
estimation variance. In this respect it is superior to 
other methods of interpolation( Journel et al 1978, 
Webster 1985, Isaaks et al 1989). 

The most common use of block kriging has been for 
the production of maps of spatially dependent proper¬ 
ties. Spatial patterns and variations can be seen clear on 
the kriging maps. Most applications of kriging in min¬ 
ing science ( Journel et al 1978, Isaaks et al 1989), soil 
science (Webster 1985, Trangmar et al 1985), ecology 
(Robertson et al 1987. 1988, Rossi et al 1992) used 
kriging maps to show spatial patterns in research re¬ 
gions or sampled locations. Block kriging is an effec¬ 
tiveness methods of interpolation for spatial data(Isaaks 
et al 1989). 

Study Site and Data 

The study site was located the old growth forests of 
Korean pine , at the Liangshui natural reserve of Xaox- 
ing'an Mountains, El28° 48’30" to El28° 55’50", N47° 
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7'39" to N47° 14'22", in Northeastern China. It is a 
typical forest vegetation and landscape in this re- 
gion(Wang et al 19%, Li et al 1993). The climate is 
terrestrial temperate and monsoon climate, average an¬ 
nual temperature is -0.3°C, and temperature range from 
-43°C to 38°C, and annual accumulated tempera- 
ture(>5°C) is about 2000°C. Annual precipitation is 
676 mm, relative humidity is 0.78, and evaporation is 
805 mm. The range of frost free season is from 100 to 
120 days. The typical soil under forests is dark-brown 
soil. 

The total area of Liangshui natural reserve is about 
6394 ha, most stands are mixture between conifer(e.g., 
pine (Pinus koraiensis), sprue e(Picea mandshurica, 
Picea koriensis ), and fir (Abies nephrolepis)) and 
hardwoods(e.g., ash, birch, oak, linden, corktree, and 
walnuts). There are twelve forest types which were dis¬ 
tributed in landscape(Li et al 1993), its age range dis¬ 
tribute from 200 to 300 years, and over 400 years for 
the dominant species of Korean pine. 

The area size for this study is 144 hm 2 , selected from 
airphoto which was taken in 1983, located in undis¬ 
turbed old growth stand. For simulating landscape pat¬ 
tern, we established 12 transects from south to north in 
1995, and setup 12 plots in each transects. The total 
numbers of plots is 144 and each plot size is 30x30 m 
in sampled region. The sampling interval between tran¬ 
sects, and between plots are 100 m(Fig.l). In each plots, 
we measured and named each trees, height, DBH, base 
area, age, cover rate, densities, soil, and plants. We also 
measured and defined the elevations, aspects, slopes, 
and micro-topography. 
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order to define and name the stand types in each plots. 
Geostatistical analysis, calculating semivariograms of 
forest types from field data, fitting various models to 
semivariograms, and kriging, were performed using 
GS+(Gamma Design 1995) software. The landscape 
heterogeneity were first analyzed by the parameters of 
semivariogram, e.g., range, sill, nugget, and fractal di¬ 
mension. The four parameters are useful tools for quan¬ 
tification of spatial heterogeneity at different landscape 
levels. Specifically, nugget variance to sill ratio reflects 
the proportion of nugget variance in total vari¬ 
ance) Robertson 1988, Li et al 1995). Where this pro¬ 
portion approaches 0, spatial dependence is very high, 
there must be spatial patterns in sampled region. Where 
this proportion approaches 1, apparent spatial depend¬ 
ence is very low, indicating that spatial pattern occurs 
mainly at scales smaller than the average distance in the 
first lag interval, in our case 100 in. 

The semivariogram of forest types, fitting by the 
models of theoretical semivariograms, can be consid¬ 
ered as landscape model for describing variations of 
forests in space in sampled region. We use ordinary 
block kriging(lsaaks et al 1989) with a block size of 20 
m across the field to interpolate unsampled locations 
and produce kriging maps. 

Results and Discussion 

Results 

There are five forest types, Spruce-fir forests) 18%), 
Spruce-Korean pine forests(48%). Libbed birch-Korean 
pine forests(9%). Amur linden-Korean pine for- 
ests(22%), Mongolian oak-Korean pine forests(3%), 
which were distributed in 144 plots and indicated dif¬ 
ferent forest ecosystems, or forest patches. 



h(m) 


Fig. 1. Diagram of plot locations sampled in the old growth 
forest of Korean pine 

In datasheet of 144 plots from fields, we calculated 
base area of each species and stand compositions, in 


Fig. 2. Semivariogram of forest types at landscape levels. The 
fitting model is linear with sill 

The variable of forest types was examined by 
semivariogram analysis. Values of the parameters of 
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semivariogram and fitting model are given in Table I. 
and isotropic semivariogram (Fig.2). four perpendicular 
directions used in anisotropic .semivariograms(Tig.3) 
were constructed. Kriging map. by using theory 
niodelf linear with sill in Table I) of semivariogram. for 
simulation landscape pattern of Korean pine forests is 
illustrated in fig.4. 


components, one can gain insight into spatial heteroge¬ 
neity at different scales. The forest types at small scales, 
independence less than 100 m. were operated by inter¬ 
processes (e.g., competition, regeneration inside stands) 
or sampling error( ITangmar et al 1085. Webster 1985). 
and the scale-dependence between 100 metros to range 
(718 m) may be caused by landscape processes. 


Discussion 

There was a high degrees of landscape heterogeneitv 
observed in old growth forests of Korean pine. 
Semivariogram for forest types show increasing!v 
strong autocorrelation a! distances of less than 718 
m(l ig.2) in landscape level. This means that forest 
patches in this landscape may he about 720 m. Beyond 
that distance the semivariogram was essentially Hat. 
indicating the region where forest types among each 
may he independence in space, the nugget vari¬ 
ance! C„-0.6200) to sill s C 1.2000) ratio show 
that the spatial heterogeneity of random is about 52%. 
and the degree of spatial dependence was moderately at 
48%. By examining the balance between these two 


• able I. Parameters of semiv ariogram and fitting model for 
forest types. 
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t ig. 4 Semis arigrams of forest tv pe in (1.45.')(! and 135 directions at landscape lev els 


I he characteristics of semivariograms of forest types 
(fig.3) in four directions show that there are a stronger 
anisotropic structure at different directions beyond the 
distances of 400 m. In direction of 0 =0. an abrupt 
change in slope of a semiv raiogram indicates that the 
variations of forest ty pes may be controlled by several 


important processes operating at different scales, while 
the semivariogram in direction of 0 -90. a flat curve 
indicates that a few processes may be equally important 
at all scales ( I rangmar et a! 1985. Webster 1985). 
Comparing fractal value (1)1.8500) with anisotropic 
structure, the landscape heterogeneity could be com- 
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plexity. Soil condition (e.g.. soil types and moisture) 
and microtopographv(e.g.. low or high site) are an im¬ 
portant factors which has effect on forest types. This is 
consistent with our Held experience in sampled region. 
Thus, the differences of anisotropic semivariograms in 
the four directions)Fig.3) should represent real differ¬ 
ences in the landscape. 




Fig. 4. A simulation of landscape pattern of old growth forest 
of Korean pine by block kirging. 

I Sprucc-tir forests. 2 Spnnv-koican pine forests. .> I ihhed hireh-kotean 
pine forests. 4 Amur liiulen-koienn pine forests. 5 Mongolian oak-korcan 
pine forests 

Landscape, or spatial, heterogeneity lead up to land¬ 
scape pattern) Forman et al 1986. Turner 1080). which 
is a scale-dependent pattern due to spatial heterogeneity 
is the function of scales. Kriging map(f'ig.4) show the 
landscape patterns of old growth forests of Korean pine 
in sampled region. From south to north, the variation of 
topography , or gradients, is about 40-60 m. Soil mois¬ 
ture is from wet in low site to dry in hill or top site. The 
lores! types or ecosystems, from spruce-fir forests in 
low-lying land to sprucc-Korean pine forests, libbed 
hireh-Korean pine forests, amur linden-Korean pine 
forests in hillside, and Mongolian oak-Korean pine for¬ 
ests in top site, were changed in series. The shape and 
size of pattern) I ig.4) show the characteristics of each 
forest types mosaic in old growth forest landscape, 
which indicate higher spatial heterogeneity. This forest 
mosaic or patches result from complex interactions be¬ 
tween pinsical(e.g.. soil, topography) and biologi¬ 
cally.g.. competition, succession) forces) Turner 1989). 


it is a unique phenomenon that emerges at the land¬ 
scape level. 

Fig. 5 illustrated the real spatial distributions of forest 
types, which were typed from datasheet of 144 plots 
and are very roughly using 100 m sampling interval. In 
kriging. there are about 3456 interpolated values with 
20 m interval. Lite field data and interpolated values 
can be used for mapping spatial pattern through kriging 
mapping. It would have some error, but we should 
know that kriging is a technique of making optimal, 
unbiased estimates of regionalized variables at unsam¬ 
pled locations using the structural properties of the 
semivaiogram and the initial set of data val¬ 
ues) I rangmar et al 1985). a useful feature of kriging is 
that an error term is calculated for each estimated value, 
providing a measure of the reliability of the interpola¬ 
tion. Comparing hotli Fig.4 and Fig.5. kriging tech¬ 
nique is an effectiveness method for simulating distri¬ 
butions of forest ty pes and to patterning in landscape. 
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Fig. 5. Spatial distributions of forests types at landscape level 
in sampled region 
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Conclusion 

Landscape ecology emphasizes broad spatial scales 
and ecological effects of the spatial patterning of eco¬ 
systems. Spatial pattern, which was operated by spatial 
heterogeneity, is the indicators of landscape and is an 
important in landscape ecological studies. Old growth 
forests of Korean pine are high degrees of spatial het¬ 
erogeneity at different scales and directions in land¬ 
scape level. The distribution of forest types in research 
region were scale-dependence and patterned in space. 
This spatial pattern which were mosiaced by forest 
types result from complex interactions between physi¬ 
cal (e.g., soil properties, topography) and biologi¬ 
cal^.g., competition, succession) forces and processes 
in landscape level. Kriging is an effective technique to 
simulate landscape patterns of old growth forests of 
Korean pine. 
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